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摘 要 在 实际 运行 过 程 中 ， 覃 式 太阳 能 集 热 器 的 聚 光 过 程 可 能 受到 太阳 光 入 射 角 、 安 装 误差 、 镜 面 误差 和 跟踪 误差 等 各 种 
实际 因素 的 影响 。 本 文 用 空间 坐标 变换 的 方法 结合 蒙特 卡 罗 光 线 追 踪 法 ， 分 别 计算 出 了 各 种 实际 因素 单独 存在 时 和 多 种 实际 因 
素 耦 合 时 的 模式 集 热 器 热流 分 布 。 结 果 显 示 ， 安 装 误差 和 跟踪 误差 使 聚焦 热流 的 不 均匀 性 增 大 ， 镜 面 误差 使 聚焦 热流 的 不 均匀 
性 减 小 ， 计 算出 的 实际 因素 耦合 热流 分 布 ， 可 以 为 槽 式 集 热 器 光学 性 能 和 热 性 能 的 分 析 提 供 依据 。 


关键 词 槽 式 集 热 器 ， 聚 光 过 程 ， 热 流 分 布 ， 光 学 误差 
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Abstract The concentration process of a parabolic trough collector can be affected by incidence angle, installation 
error, reflector slope error and tracking error in operation. The coordinate transformation method and the Monte 
Carlo Ray Trace method are integrated to compute heat flux distributions under individual and coupled optical errors 
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0 引言 1.1 热流 分 布 模拟 方法 
模式 太阳 能 集 热 器 (PTC) 是 太阳 能 热 发 电 的 Sg a te es 


种 重要 方式 ， 在 对 PTC 进行 研究 时 ， 集 热管 外 聚焦 ” 画 出 玻璃 管 。 反 射 镜 张 口 宽度 、 焦 距 和 集 热 需 长 度 
热流 的 计算 常常 是 第 一 步 。 何 雅 玲 中 使 用 蒙特 卡 罗 光 分 列 用 wf 和 工 表示 ,坐标 原 上 扩建 并 在 吸 热 管 的 进 
线 追 踪 法 (MCRT) 计算 了 理想 情况 下 PTC 的 聚焦 热 。 口 截 面 国 心 。 

流 ， 证 明了 MCRT 方法 的 实用 性 ，MCRT 方法 得 到 NE 
了 广泛 的 应 用 P21。 与 理想 情况 相 比 ，PTC 的 聚 光 过 ”假如 入 射 光线 万 向 向 量 为 M ， 反 射 面 的 法 向 量 为 


程 受到 各 种 实际 因素 包括 太阳 光 入 射 角 、 安 装 误差 、 。 内 ， 那 么 反射 光线 方向 向 量 思 的 计算 公式 为 ; 

跟踪 误差 和 镜面 误差 的 影响 , 赵 东 明 等 "采用 空间 华 六 入 办 这 a) 
te te A 各 种 实际 因素 实际 上 是 对 集 热管 位 置 、 入 射 光 
人 线 或 镜面 法 线 方向 产生 影响 ， 在 MCRT 的 算法 中 ， 


寸 镜 面 误 弛 定 误 差 进 详尽 的 敏感 分 Ns > ~» » 
对 镜面 误差 和 跟踪 误差 进行 了 详尽 的 敏感 性 分 析 ， 这 三 者 间 互 不 干涉 。 为 了 研究 各 种 因素 的 影响 ， 只 
均 发 现 当 这 些 实际 因素 存在 时 ， 聚 焦 热 流 将 发 生 很 需要 改变 集 热管 的 位 置 ， 或 用 空间 坐标 变换 的 方法 
后 西 ee 而 女 厅 和 fy 上 业 且 ， 疏 用 工 小 万 + 

大 的 改变 ， 因 此 需要 对 各 种 实际 因素 存在 时 PTC 的 。 计算 新 的 商 或 方 即 可 。 
集 热管 外 聚焦 热流 进行 研究 。 

目前 的 文献 常常 只 计算 出 只 有 一 种 实际 因素 存 
在 时 的 PTC 聚焦 热流 ， 而 实际 应 用 中 ，PTC 的 聚 光 
过 程 可 能 会 同时 受到 多 种 实际 因素 的 影响 。 本 文 先 
计算 各 种 实际 因素 单单 独 存在 时 PTC 的 聚焦 热流 3} 
进行 分 析 ， 然 后 在 MCRT 方法 的 基础 上 使 用 空间 坐 
标 转换 ， 计 算出 多 种 实际 因素 耦合 时 的 聚焦 热流 。 
1 模型 描述 (a) 截面 示意 图 
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(b) 三 维 视图 

图 1 覃 式 集 热 器 示意 图 

Fig. 1. Schematic of a PTC (a) cross sectional view (b) 3D 
View 


1.2 安装 误差 
当 集 热管 的 安装 误差 使 级 热管 的 位 置 发 生平 移 
时 ， 假 设 管子 圆心 线 在 y 和 z 方向 新 的 坐标 分 别 为 
yj 和 zy， 则 吸 热管 外 壁 方程 变 为 : 
(y-») +(z-a) =7 O) 
r 是 吸 热管 的 外 壁 半径 。 集 热管 安装 误差 或 者 支 
撑 架 的 形变 也 可 能 会 造成 集 热管 的 位 置 与 焦点 线 不 
平行 。 当 PTC 采用 双 轴 跟踪 系统 时 ， 光 子 被 吸 热管 
接收 位 置 的 x 坐标 等 于 镜面 上 入 射 点 的 x 坐标 《用 
2 表示 )， 所 以 只 要 知道 吸 热管 中 心 线 的 坐标 方程 ， 
光子 到 达 吸 热管 的 位 置 仍然 可 以 求 出 来 。 
简单 起 见 ， 假 设 吸 热管 的 进口 截面 圆心 坐标 为 
(0.0.0), 出 口 埠 面 圆心 坐标 为 (yy 2， 那么 在 反射 】 
线 所 在 的 垂直 于 焦点 线 的 平面 内 ， 吸 热管 外 壁 的 4 
标 方程 为 : 


2 2 
人 ~ 
人 + -至 四 


将 吸 热管 新 的 坐标 方程 (2) 或 (3) 代 入 MCRT 算法 
中 ， 保 持 其 他 部 分 不 变 ， 即 可 求 得 有 集 热管 安装 误 
差 的 管 外 不 均匀 热流 分 布 。 

对 于 单 轴 跟 踪 系 统 ， 由 于 入 射 角 的 存在 ， 反 射 
光线 和 集 热 管 外 壁 的 交点 不 易 确 定 ， 情 况 更 为 复杂 ， 
本 文 不 作 研 究 。 
1. 3 太阳 光 入 射 角 

PTC 对 太阳 位 置 的 跟踪 有 单 轴 跟踪 和 双 轴 跟踪 


卜 3 


两 种 ， 只 有 双 轴 跟踪 系统 才能 使 太阳 光线 保持 直射 
入 射 ， 但 是 由 于 双 轴 跟踪 系统 结构 复杂 ， 工 程 实际 
中 PTC 一 般 采 取 单 轴 跟 踪 系 统 ， 这 种 情况 下 入 射 太 


阳 入 射 角 往 往 不 为 0"， 如 图 2 所 示 。 
太阳 光线 4B 以 入 射 角 9 照射 到 反射 镜 边缘 上 的 
点 了 ,然后 被 反射 到 吸 热管 上 的 点 C， 距 离 CD 即 为 
经 过 B 点 的 太阳 光线 不 能 照射 到 的 长 度 。 文献 求 出 

CD 的 最 大 值 和 最 小 值 分 别 为 

要 十 次] tan O (4) 
n=f tan0 (5S) 
Lmax 是 契 PIC a 册 损 失效 应 色 By 影响 到 的 最 大 距 
离 ， 即 与 末端 D 距离 大 于 Zaax 的 区 域 总 的 聚焦 热流 
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不 变 ， 而 与 末端 D 距离 小 于 Lmin 的 管 段 则 不 能 接收 
到 反射 后 的 太阳 光 。 
2 六 ”太阳 


图 2 太阳 光 斜 射 入 射 图 
Fig. 2. Diagram of an incident ray on the sun-facing end of the 
reflector 


图 3(a@) 中 光线 1 是 MCRT 法 计算 出 来 的 直射 入 
射 光线 ， 光 线 2 是 入 射 角 为 9 的 斜 射 入 射 光线 。 图 
3(b) 中 的 坐标 系 O-xyg 是 坐标 系 O-xyz 绕 y 轴 旋 转角 
度 9 后 得 到 的 ， 光 线 1 在 此 坐标 系 中 的 坐标 就 是 光 
线 2 在 坐标 系 0-xtyz 中 的 坐标 ， 用 M 在 表示 ， 则 | 
M'’=A:M (6) 
菱 中 4 是 从 坐标 系 O-xyz 到 坐标 系 O-xyz 的 空 
间 坐 标 变换 矩阵 ; 


coSsO_ 0 SnO 
A=| 0 1 0 (7) 

—sin0 0 cosO 
将 得 出 的 M 代 入 公式 (D 中 ， 保 持 其 他 部 分 不 
变 , 即 可 求 得 太阳 光 入 射 角 不 等 于 0" 时 的 热流 分 布 。 


y (y') 


人) 空间 坐标 变换 
3 单 轴 系 统 太 阳光 入 射 角 坐 标 计 算 
Fig. 3. The calculation of the incident direction of a single-axis 
PTC (a) different incident conditions where ray 1 is the vertical 
incident ray and ray 2 is the oblique incident ray (b) the 
coordinate transformation for the incidence angle 0 


1.4 J 

跟踪 误差 是 由 于 跟踪 不 准确 导致 的 集 热 器 位 置 
的 偏差 ， 如 图 4(a) 所 示 。 由 于 跟 踩 误差 的 存在 ， 集 热 
器 1 相对 于 正确 的 位 置 ( 集 热 器 2) 逆 时 针 旋转 了 
角度 。 如 图 4(b) 所 示 ， 空 间 坐 标 变换 矩阵 等 于 : 


1 0 0 
于 cosp Sn (8) 
0 -sinp cosp 
则 修正 后 的 入 射 光 线 坐 标 M' 为 : 
M'’=B:M (9) 
aly 保持 其 他 部 分 不 


变 ， 即 可 求 得 有 跟踪 误差 的 热 
海 


\ 流 分 布 。 


M (1) 


(b) 空间 坐标 变换 
图 4 跟踪 误差 及 其 坐标 变换 
Fig. 4 Tracking error and its coordinate transformation (a) 
diagram of tracking error (b) the coordinate transformation for 
the tracking error fp 


1.5 镜面 误差 

镜面 误差 考虑 的 是 反射 面 微 观 的 缺陷 或 宏观 的 
变形 等 原因 所 造成 的 反射 误差 。 镜 面 误差 包含 的 种 
类 很 多 , 把 所 有 的 镜面 误差 锥 加 起 来 后 可 以 用 实际 
的 与 理想 的 镜面 法 线 方向 的 偏差 来 代 蔡 趾 , 认为 此 偏 
差 角度 在 x 方向 和 y 方 向 分 解 后 得 到 的 两 个 分 量 分 别 


符合 相同 的 高 斯 分 布 四 ， 镜面 实际 法 线 坐标 向 量 N' 
可 以 由 理想 的 镜面 法 线 方向 N 经 过 两 次 坐标 变换 得 
到 : 
N'=B.(A:N) (10) 
其 中 4 和 B 用 公式 (7) 和 (8) 来 计算 ， 公式 中 的 角 
度 9 和 的 概率 密度 函数 相同 : 
es | (11) 
J ge 
1 xXx? 
ee | (12) 


o 是 标准 方差 , 也 是 镜面 误差 角度 。 将 得 出 的 N， 
代入 公式 (1) 中 ， 保 持 其 他 部 分 不 变 ， 即 可 求 得 有 镜 
面 误差 的 热流 分 布 。 

1.6 实际 因素 耦合 


V 会 人 乍 期 于 | 
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对 于 安装 误差 ， 与 理想 模型 相 比 ， 吸 热管 坐标 
发 生 了 变化 ， 用 式 (2) 或 式 (3) 来 考虑 其 影响 即 可 。 对 
于 跟踪 误差 、 镜 面 误差 、 太 阳光 入 射 角 等 实际 因素 ， 
可 以 用 以 下 公式 计算 反射 光线 方向 : 

P=-2(N'-M M’')N'+M (13) 

其 中 NW' 和 闻 分别 是 考虑 了 跟踪 误差 、 镜面 误差 

和 太阳 光 入 射 角 等 因素 的 实际 反射 镜 法 线 和 实际 入 

射 光线 坐标 向 量 。 

需要 注意 的 是 ， 由 于 和 矩阵 乘法 不 满足 交换 律 ， 

所 以 坐标 变换 的 顺序 不 是 随意 的 。 例 如 ， 当 同时 考 

虑 入 射 角 影 响 和 跟踪 误差 时 ， 实 际 入 射 光线 坐标 
为 : 


i 


Be B.(A: M) (14) 
采用 M'= 8.( A: 呈 或 M'=A.(B-MM) 得 到 的 
结果 并 不 一 样 。 ee 


光线 ， 甜 阵 考虑 跟踪 误差 ， 考 虑 到 实际 过 程 ， 应 
该 采取 M'= B.(4.M) 来 计算 跟踪 误差 存在 时 单 轴 
pe a 


结果 与 讨论 

下 面 利用 空 de 
阳光 入 射 角度 等 实际 因素 分 别 进行 模拟 和 分 析 ， 
后 将 多 种 实际 因素 耦合 在 一 起 计算 。 因为 管 外 热流 
分 布 的 相对 大 小 与 直射 辐射 强度 没有 关系 ， 所 以 分 
析 结 果 时 管 外 热流 取 的 是 当地 聚焦 比 〈LCR )。 本 文 
模拟 时 采用 LS-3 型 模式 太 阳 能 集 热 器 作为 物理 模 
型 ， 名 尺 3 车 参数 见 表 1。 关 于 其 他 型 号 集 热 器 的 


结论 可 以 通 过 用 相同 的 方法 分 析 得 到 。 
表 1 LS-3 模式 集 热 器 数据 

Table1 Parameters of the LS-3 PTC 
参数 数值 
吸 热管 半径 mm 0.035 
焦距 /m 1.71 
反射 镜 张口 宽度 /m 5.76 
玻璃 管 透射 率 0.95 
吸 热管 吸收 率 0.96 
反射 镜 反 射 率 0.93 


2.1 安装 误差 的 影响 

y 方 向 和 z 方 向 安装 误差 的 定义 分 别 为 :Y=yj/w 
和 2Z=zy/f。 图 5 中 显示 的 是 平移 安装 误差 下 的 热流 分 
布 。 集 热管 安装 误差 达到 毫米 数量 级 时 ， 管 外 热流 
形状 就 已 发 生 明 显 改变 。 图 6 和 图 7 比较 了 两 方向 
安装 误差 都 存在 时 y 或 z 方 向 安装 误差 对 管 外 热流 分 
布 的 影响 。 图 6 保持 Z 不 变 改变 7 了 的 大 小 ， 热 流 分 
布 形状 的 改变 非常 大 ， 图 7 保持 了 不 变 改变 Z 的 大 
小 ， 热 流 分 布 形状 的 改变 相对 较 小 。 因 此 , 与 z 方 向 
安装 误差 相 比 ，y 方向 安装 误差 的 影响 大 。 所 以 ,在 
集 热管 的 安装 过 程 中 要 格外 注意 减少 在 y 方向 的 安 
装 误差 。 


(b) 2Z=0.5% 
图 62Z 相同 时 不 同 了 对 聚焦 热流 的 影响 
Fig. 6 Flux distributions for different Y when Z is constant (a) 
2Z=0.2% (b) 2Z=0.5% 
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(a) z 方 向 安装 误差 
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图 5 平行 安装 误差 下 的 热流 分 布 0 一 一 上 
Fig. $5 Flux distributions for installation errors (a) z direction 0 60 120 180 240 300 360 
installation errors (b) y direction installation errors 周 向 角度 
140 人 (b) 7¥=0.5% 
BE 7 了 相同 时 不 同 2Z 的 影响 
120 | | 03% Fig. 7 Flux distributions for different Z when Yis constant (a) 
1 04% 
100| | [os Y=0.2% (b) Y=0.5% 
本 
0 | 图 8 显示 的 是 不 平行 安装 对 热流 分 布 的 影响 。 
D ee i Ee 、 BA 
60| A | ] 集 热管 进口 截面 圆心 坐标 为 (0,0,0)， 出 口 截面 圆心 坐 
40| A | ll | 标 为 (4,y1,z/)， 其 中 六 和 zj 分 别 对 应 天 0.1% 和 
岂 | > 、 Ss 
下 Z-0.1%。 昌 然 图 中 热流 分 布 在 轴线 上 是 不 均匀 的 ， 
20 上 嘱 l a > ee ,EA 、A- 
有 记 但 可 以 看 成 无 数 个 平移 安装 误差 的 热流 分 布 连续 排 
00 60 120 180 240 300 360 布 而 成 。 
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周 向 角度 Fig. 8 Flux distributions in case of Y=0.1% and 2Z=0.1% for 


outlet of the absorber tube 


2.2 入 射 角 的 影响 

图 9(a) 显 示 的 是 直射 辐射 强度 相同 时 双 轴 跟踪 
和 单 轴 跟 踪 系 统 热 流 沿 管 长 的 分 布 ， 纵 轴 为 截面 吸 
收 的 总 能 量 。 对 于 单 轴 跟 踪 的 系统 ,图 中 画 出 了 Lnax 
和 Lwin 的 位 置 ， 可 以 看 到 MCRT 方法 模拟 的 结果 和 
公式 计算 结果 吻合 得 很 好 ， 证 明了 MCRT 方法 对 太 
阳光 和 斜 射 入 射 时 热流 分 布 模拟 的 准确 性 。 在 不 受 末 
集 热 管 每 个 截面 接收 的 总 
的 聚焦 热流 等 于 太阳 光 以 相等 的 直射 辐射 强度 直射 
Www 
是 不 变 的 。 A A 


热管 接收 到 的 聚焦 能 量 从 零 开始 连续 地 增 大 ， 最 终 
成 为 周 向 不 均匀 、 多 向上 下 要 的 全 
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(b) 聚焦 热流 三 维 视图 
图 9 太阳 光 30° 入 射 角 热 流 分 布 
Fig. 9 Flux distributions at 30° incidence angle (a) axial 
distribution (b) 3D distribution 


2. 3 镜面 误差 和 跟踪 误差 影响 

图 10 显示 的 是 镜面 误差 对 聚 光 过 程 的 影响 。 从 
图 10(a) 可 以 看 出 ， 周 向 热流 分 布 仍然 是 对 称 的 ， 镜 
面 误差 的 作用 是 使 聚焦 热流 由 不 均匀 趋向 均匀 ， 这 
种 作用 有 利于 减 小 管 壁 温差 ， 且 镜面 误差 越 大越 明 
显 。 图 10(b) 显 示 的 是 镜面 误差 角度 对 未 受 末 端 损 失 
影响 区 域 的 拦截 系数 的 影响 ， 当 镜面 误差 大 于 0.12° 
时 ， PTC 的 拦截 系数 小 于 1， 应 该 避免 这 种 情况 。 
图 11 显示 的 是 跟踪 误差 对 聚 光 过 程 的 影响 。 从 
图 11(a) 可 以 看 出 ， 跟踪 误差 使 太阳 光 1 离 垂 直 于 反 
射 镜 张 口 平面 的 方向 入 射 ， 增 大 了 热流 的 不 均匀 性 。 
图 11(b) 显 示 的 是 跟踪 误差 角度 对 未 受 末端 损失 影响 
区 域 的 拦截 系数 的 影响 ， 当 跟踪 误差 大 于 0.6° 时 , 拦 


截 系数 小 于 1。 
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图 10 镜面 误差 对 聚 光 过 程 的 影响 

Fig. 10 Effects of reflector Slope error on concentration process 
(a) circumferential heat flux distributions (b) variation of the 

intercept factor 
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Fig. 11 Effects of tracking error on concentration process (a) 


circumferential heat flux distributions (b) 
intercept factor 


2.4 实际 因素 耦合 


呈 


太 
为 
并 
可 


variation of the 


LCR 
60.00 


S0.00 
40.00 
30.00 
20.00 
10.00 


0.000 


LCR 


20 上 


0 60 120 180 240 300 360 


周 向 角度 
(b) 周 向 热流 分 布 


12 太阳 光 入 射 角 30"*，7=0.1%，2=0.1%， 跟 踪 误 差 0.2°， 
镜面 误差 0.1° 时 的 热流 分 布 

Fig. 12 Flux distributions in case of 30° incidence angle, 

¥=0.1%, 2Z=0.1%, 0.2° tracking error and 0.1° reflector Slope 

error (a) 3D distribution (b) circumferential distribution 


图 12 显示 了 多 种 实际 因素 耦合 


时 的 热流 分 布 ， 


阳光 入 射 角 为 30°*，7 了 =0.1%，2Z=0.1%， 跟 踪 误 差 


0.2*， 镜 面 误差 为 0.1°。 
损失 影响 区 域 的 周 向 热流 。 由 


图 12(b) 显 示 的 是 未 受 末 
于 各 种 实际 误差 角 


的 大 小 不 一 ， 所 以 管 外 聚焦 热流 的 形状 变化 多 样 ， 
以 计算 出 集 热 管 外 实际 因素 耦合 的 热流 分 布 ， 从 


而 为 管内 三 维 温度 场 的 计算 提供 依据 。 
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结 论 


本 文 分 别 分 析 了 集 热 管 安装 误差 、 


镜 


响 
得 
时 


流 
向 
注 
能 


集 
下 


bE 量 沿 轴 


太阳 光 入 射 角 、 


误差 和 跟踪 误差 等 实际 因素 对 PTC 热流 分 布 的 影 


， 计 算出 了 多 种 实际 因素 耦合 时 的 PTC 热流 分 布 ， 


出 了 以 下 结论 ; 


1) 集 热 管 安 装 误 差 达到 毫米 数量 级 时 ， 聚 焦 热 
形状 融 已 发 生 明 显 改 变 ，y 方向 的 安装 误差 比 z 方 
对 热流 分 布 影响 更 大 ， 集 热管 安装 过 程 中 要 格外 


意 减 少 在 > 方向 的 安装 误差 ; 
2) 在 末端 损失 影响 区 域 ， 吸 热 
向 从 零 开 始 连 续 地 增 大 ， 


管 接收 到 的 聚焦 
在 此 区 域 以 外 ， 


热管 每 个 截面 总 的 聚焦 热流 等 于 太阳 光 以 相等 的 
直射 罚 射 强度 直射 入 射 时 总 的 聚焦 热流 ， 且 在 轴 向 


可 以 认为 不 变 ; 


均匀 ， 这 种 作 


用 误差 的 作用 是 使 聚焦 热流 由 不 均匀 趋向 
有 利于 减 小 管 壁 温差 ,镜面 误差 越 


3) 镜 


大 越 明显 ， 而 跟踪 误差 增 大 聚焦 热流 的 不 均匀 性 ， 


= 
忒 0 


响 集 热管 的 安全 性 ， 当 镜面 误差 大 于 0.12? 或 跟踪 


误差 大 于 0.6* 时 ，PTC 的 拦截 系数 小 于 1; 
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(4) 由 于 各 种 实际 误差 角度 的 大 小 不 一 ， 所 以 


车 外 聚焦 热流 的 形状 变化 多 样 ， 可 以 计算 出 集 热 管 


下 


外 实际 因素 耦合 的 聚焦 热流 ， 从 而 为 光学 性 能 和 热 
生 能 的 分 析 提供 依据 。 
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